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Abstract 

We investigate the T(3)-gauge theory of static dislocations in continuous solids. 
We use the most general linear constitutive relations bilinear in the elastic distor- 
tion tensor and dislocation density tensor for the force and pseudomoment stresses 
of an isotropic solid. The constitutive relations contain six material parameters. In 
this theory both the force and pseudomoment stresses are asymmetric. The theory 
possesses four characteristic lengths t\, £2, £3 and £4 which are given explicitely. 
We first derive the three-dimensional Green tensor of the master equation for the 
force stresses in the translational gauge theory of dislocations. We then investigate 
the situation of generalized plane strain (anti-plane strain and plane strain). Using 
the stress function method, we find modified stress functions for screw and edge 
dislocations. The solution of the screw dislocation is given in terms of one indepen- 
dent length £\ = £4. For the problem of an edge dislocation, only two characteristic 
lengths £2 and £3 arise with one of them being the same £2 = £\ as for the screw 
dislocation. Thus, this theory possesses only two independent lengths for gener- 
alized plane strain. If the two lengths £2 and £3 of an edge dislocation are equal, 
we obtain an edge dislocation which is the gauge theoretical version of a modified 
Volterra edge dislocation. In the case of symmetric stresses we recover well known 
results obtained earlier. 

Keywords: Gauge theory; dislocations; defects; incompatible elasticity, torsion, 
Green tensor. 
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1 Introduction 



1.1 Reviewing remarks 

Gauge theories are very successful in high energy physics to understand and to explain 
fundamental interactions (see, e.g., Important solutions in gauge theories are 

topological strings like the so-called Abrikosov-Nielson-Olesen string 0,13]. It is also 
known that the theories of gravity and generalized gravity can be understood as gauge 
theories of the Poincare group and the affine group, respectively 0, @]. For example, 
the teleparallel formulation of Einstein's general relativity can be understood as a gauge 
theory of the four-dimensional translation group (see, e.g., Hehl [7|). Thus, the torsion 
tensor is not just a tensor, but rather a very special tensor that is related to the translation 
group . In such theories mathematical solutions for cosmic strings and other topological 
defects in space-time are found [§]. But it is hard to observe them by experiment. 

An application of the three-dimensional translational gauge theory is the gauge theo- 
retical formulation of dislocations. The torsion tensor plays the role as the translational 
field strength. In such a theory dislocations appear quite naturally as topological defects 
(strings). Dislocations play an important role in solid state physics |9[]. It is easy to 
investigate them experimentally. Therefore, the gauge theory of dislocations is the best 
example to investigate the translational gauge theory as a physical field theory. Pioneering 
work in the field of the translational gauge theory of defects in solids was made by Kadic 
and Edelen Edelen and Lagoudas [12|. It is well known that dislocations are the 



fundamental carrier of plasticity. The hope of the gauge theory of dislocations is to un- 
derstand plasticity as a 'fundamental' interaction in solids. Materials display strong size 
effects when the characteristic length scales are of the order of micrometers down to a few 
nanometers. With the interest in miniaturization of devices, the length scales associated 
with nanodevices are very small and classical theories of elasticity and plasticity are not 
applicable. Classical elasticity and plasticity theories cannot explain the size dependence 
because their constitutive laws possess no internal material lengths. In metals, the phys- 
ical origin of size effects can be related to dislocations. Gauge theory of dislocations may 
be used to explain size effects and to model mechanical properties of miniaturized devices. 
In the gauge theory of dislocations the dislocation density tensor which is the curl of the 
elastic distortion or of the negative plastic distortion multiplied by material parameters 
which define internal material lengths enter the constitutive relations. The gauge theory 
of dislocations shall be a step in the direction of a microscopic theory of plasticity that 
incorporates interaction of structural defects. The dislocation should play a similar role 
as a gauge boson like the photon as a gauge boson in Maxwell's theory. Nevertheless, 



the theory of Kadic and Edelen Edelen and Lagoudas [12[ has some lacks. They 
used very special constitutive relations. For example, the force stress tensor is symmetric. 
But in the continuum theory of dislocations it is known that the force stress tensor has 



to be asymmetric at least in the dislocation core region [13l . Il4| and that moment stress 



is the specific response to dislocations [15] . Such a moment stress cannot be calculated 
from elastic moduli only. A realistic constitutive law between the moment stress and the 
dislocation density tensor has to be used to built a physical dislocation theory. 

Therefore, such a theory is not the most general isotropic dislocation gauge theory. 
Their theory is just able to give an acceptable solution of a screw dislocation after some ad 



hoc assumptions [ly, [17( . Moreover, it is not possible to give a correct solution of an edge 



dislocation. Especially, the stress component a zz is incorrect and the condition of plane 
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strain is not satisfied 18j. Such a dislocation gauge theory with symmetric force stresses 
possesses only one additional material parameter and two length scales can be defined to 
describe size-effects. In order to explain bending and torsion experiments one needs at 
least two length scales. One question that arises is how many characteristic length scales 
should appear in an improved gauge theory of dislocations? Recently, Sharma and Ganti 
19( | have combined Edelen's gauge model with a simple strain gradient elasticity. These 
authors applied it to a screw dislocation. 

These drawbacks of the Edelen model have been a motivation for further investiga- 



tions 
Lazar 
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The so-called Einstein choice has been investigated by Malyshev pi| and 
Their solution of a screw dislocation looks to be acceptable for theories 
with symmetric force stresses. In the solution of an edge dislocation found by Malyshev 
18j, a strange far field emerges, which shows an oscillating behavior. In addition, the 
condition of plane strain is not fulfilled. Later, Lazar [23[ tried to find a more realistic 
gauge theoretical solution for an edge dislocation. He used another constitutive relation 
and it was necessary to introduce additional gradients of the elastic rotation which is a 
little sophisticated. Thus, up to now, the case of an edge dislocation has not been treated 
satisfactory in the gauge theoretical framework. Quite recently, Malyshev 24| has used 
the Einsteinian T(3)-gauge approach to calculate the stress tensor of a screw dislocation 
in the second order. Some other proposals are given by Katanaev and Volovich 25j, Popov 
and Kroner 



23 



Kroner 27 



We want to note that all these kinds of gauge theories of 
dislocations use symmetric force stresses and some special choices for the constitutive law 
of the moment (or hyper) stresses. 



1.2 The subject of the paper 

In this paper we want to investigate the translational gauge theory of dislocations with 
the most general linear isotropic constitutive relations quadratic in the physical state 
quantities for the force and pseudomoment stress tensors with six material coefficients. 



Such a theory was originally proposed by Lazar 21]. However, Lazar [2JJ used a simplified 
version for dislocations. Our purpose is to derive the general framework of such a theory 
and to apply it to find new analytical solutions of straight screw and edge dislocations 
in an infinite medium. The paper consists of five sections. Section [T] is the introduction. 
In section [2] we give the basics of the translational gauge theory of dislocations. We 
give the relation for positive definiteness of the elastic energy and the dislocation core 
energy. We give the relation between the pseudomoment stress tensor and the moment 
stress tensor. We derive the Euler-Lagrange equation in terms of the asymmetric stress 
tensor. This is the central equation for an isotropic material with dislocations. We set 
up the relation to other gauge theories of dislocations proposed earlier. We calculate the 
three-dimensional Green tensor for the central equation of the asymmetric force stresses. 
In this manner we are able to introduce four characteristic length scales of the gauge 
theory with asymmetric force stresses. Later we investigate dislocations in the framework 
of generalized plane strain. Sections [3] and @] present the gauge theoretical solutions of 
straight screw and edge dislocations in an infinite medium with asymmetric force stresses. 
We will calculate the elastic distortion and dislocation density (torsion) tensors and the 
force and pseudomoment stresses of screw and edge dislocations. We prove that the force 
stresses of a screw dislocation and an edge dislocation have to be asymmetric instead 
of symmetric in order to guarantee nonnegative energy. Section concludes the paper. 
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Details of dislocations in asymmetric elasticity are given in the appendices \K\ and [Bj 



2 Gauge theory of dislocations 
2.1 Foundations 

In the translational gauge theory of dislocations the displacement vector Ui is defined up 
to a gauge transformation - a local translation Ti(x) 

Ui^>Ui + n(x) (2.1) 

and the translational gauge field <f>ij transforming under the local translation in a suitable 
form 

</>ij -»• faj - Tij(x). (2.2) 

From the mathematical point of view, the gauge field <pij is the translational part of the 
generalized affine connection 2(|. In the translational gauge theory, the elastic distortion 



is defined by the translational gauge- covariant derivative according to 

Pij = Y? M i := u iJ + far ( 2 - 3 ) 

Thus, the elastic distortion is gauge independent and, therefore, a physical state quantity. 
The asymmetric tensor 0jj may be identified as the negative plastic distortion tensor [20, 



211 ] . The elastic strain energy is a function of the elastic distortion. It reads 



with the specific response 



W d = i/Vfc ( 2 - 4 ) 



9Wa (o *\ 



It is the analog of the Piola-Kirchhoff stress tensor. It is an asymmetric stress tensor. 
The most general isotropic constitutive relation has three material coefficients. Thus, the 
asymmetric force stress tensor has the form0 



a 



ij 



A Sijp kk + 2/j,p (ij) + 27/3^]. (2.6) 



Here \i and A are the Lame coefficients. The coefficient 7 is an additional material 
parameter due to the skew-symmetric part of the elastic distortion (the elastic rotation). 
Thus, 7 is the module of rotation (see also Kroner [28]). The skew-symmetric stress ofei is 



caused by the (local) elastic distortion (3^. If 7 = in Eq. (|2.6p . we obtain the Hooke law 
with <T[ij] = 0. Thus, in this case the elastic rotation (3^ is undetermined. The positive 
semi-definiteness of the elastic distortion energy W e \ > requires the restriction 

fi>0, 7>0, 2/i + 3A>0. (2.7) 



1 We are using the notations A/ij) = + Aji) and Auj] = h(Aij — Aji). 
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Another physical state quantity is the translational gauge field strength - the torsion 
tensor 

Tijk 4*ik,j 4*ij,kj T^ijk Tikj (2-8) 

and 

Tijk ftikj Pij,k- (2-9) 

It has the physical meaning of the dislocation density tensor 

c% = -tjkiTiki = £jki(3u,k- (2.10) 
They fulfill the translational Bianchi identity 

£jkiTijk,i = 0, a ij:j = 0. (2.11) 

Thus, the physical state quantities are the elastic distortion, and the dislocation 
density tensor T^ k . 

The dislocation core energy contains only the translational gauge field strength and is 
given by 

Wdisi = - TijkHijk (2.12) 
with the translational gauge field momentum tensor 

Hijk = 2-^= — , Hij k = —H ik j. (2-13) 

01 ijk 

The translational gauge field momentum tensor is a very special hyperstress tensor. 
It is the response quantity to the dislocation density tensor and has only nine independent 
components. It is convenient to perform an irreducible decomposition of the torsion, from 
which we can construct the most general linear and isotropic constitutive law between the 



dislocation density and the hyperstress in the following way [20l . |21 

3 



H l]k = J2 a i {I %k- (2.14) 



i=i 



Here a±, a,2 and 03 are nonnegative material coefficients due to the positive semi-definiteness 
if the dislocation core energy VFdisi > 0. They have the dimension of a force. The irre- 
ducible decomposition of the torsion under the rotation group 5*0(3), with the numbers 
of independent components 9 = 5 © 3 © 1, is given by 

Tijk ^ ^Tijk © ^ ^Tijk © ^ ^Tij]~. (2.15) 

The tensor, the trace and the axial tensor pieces are defined by (see also 0, [2l| ) 

(1 % k = T ijk - &>T ijk - ^T ijk (tentor) (2.16) 

{2) T ijk := - (SijTuk + 5 ik Tiji) (trator) (2.17) 

(3) T ijk : = i (T ijk + T jki + T kij ) (axitor). (2.18) 
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The hyperstress tensor may be written as 



Hijk — Ci Tij k + c 2 (Tjki + T ki j) + c 3 [SijTuk + S ik Tiji) 
with the abbreviations 



(2.19) 



ci := -(2oi + o 3 ), <fc := - (a 3 - ai) , c 3 := -(a 2 - ai). (2.20) 



The inverse relations are 



a\ — ci — c 2 > 0, a 2 = ci - c 2 + 2c 3 > 0, a 3 = 2c 2 + Ci > 0. 



(2.21) 



Since and are uncoupled from each other, the conditions of positive semi-definiteness, 
W = W e \ + Wdisi > 0, can be studied separately: W e \ > and Wdisi > 0. Thus, the char- 
acteristic constants of the material have to satisfy the conditions (12.71) and (I2.21j) . 

With the relations H im = \ej km Hijk and = €j kn ai n , we obtain from Eq. (12.141) the 
irreducible decomposition 



Hij = ai ( «(y) - - 5ija kk j + a 2 a^-j + a 3 - <%a fcfc 

trator: 



tentor: 



(2.22) 



Thus, the axitor represents a set of three perpendicular forests of screw dislocations of 
equal strength. The trator describes forests of edge dislocations: a 



-a. 



and a 



-a 



The tentor is the deviator of the dislocation density tensor and it 



represents forests of edge dislocations: a a 
dislocations: 



(i) 



) with (V +« 



0, 



(1) 



2,^ ^yxi ^xz "zi 

3 (c^xa; "I - Cl^y "I - Ct zz y n u-n n.r.r > ' ///; 

relation was given earlier by Hehl and Kroner [13[], Kroner 



a ?T and a 



a, 



OLfiq, 



a zy and screw 

'mi 3 1 1 '-xx ^yy ^zz) : 

= 0. A similar constitutive 
14 , Katanaev and Volovich 



251 ]. Popov and Kroner [26j|. Alternatively, the dislocation core energy (I2.12p may be 



rewritten in the form 



disl 



2 a ijHij- 



[2.23) 



We want to note that a similar hyperstress tensor like (I2.22p in terms of the dislocation 



density tensor and with three material coefficients was derived by Eringen and Claus 29 



111 



a special case of micromorphic elasticity. Therefore, the static gauge theory of dislocations 
is related to a micromorphic theory and not to a Cosserat theory if we identify the gauge 



field 4>ij with the micro-distortion tensor of micromorphic elasticity [30l . l3l| . In general 



a hyperstress tensor possesses 27 independent components and is the response quantity 
to the gradient of the micro-distortion tensor. The isotropic constitutive relation of such 
hyperstress tensor has 11 material parameters (see, e.g., Mindlin [30|). 

Since the tensor Hiuu has only nine independent components like a couple stress 
tensor, the following question arises: In which way is the tensor Hij k related to the couple 
stress tensor? The contortion tensor is defined in terms of the torsion tensor according 
to 



32, 33, 34 



Kij k 



" 2 \Tijk + Tjki 



Kij k Kji k . 



(2.24) 



6 



A moment stress tensor possesses nine independent components and for an isotropic ma- 
terial the constitutive relations has just three material parameters. The contortion tensor 
contains nine independent rotational degrees of freedom. Obviously, is antisymmet- 
ric in i, j, whereas T^k is so in j, k. With the identity Ki 
so-called Nye curvature tensor |35J (see also |36|, |37J, |38|, [39 



\eijiKijk, we recover the 



a 



31 2 $ij a lh 



atij 



(2.25) 



The moment stress (or spin moment stress tensor) is the response of the crystal to the 
contortion 



dW 

dKi jk '' 



[2.26) 



Thus, contortion produces moment stress and vise versa. We can rewrite the tensor 
in terms of the moment stress tensor as follows 



H, 



ijk 



[ij]k 



7~ijk Tjki 



(2.27) 
(2.28) 



Therefore, the tensor is given in terms of moment stresses similar like the torsion 
in terms of the contortion. Such a tensor H^u is called pseudomoment stress tensor or 
sometimes spin energy potential tensor 41, 42, i|, @|. In addition, we have 

1 



H 



i.i 



• 1.1 



Hji — SijHu, 



(2.29) 



where = dW/dKij . We conclude, torsion produces pseudomoment stress and vise 
versa. Thus, pseudomoment stress is the specific response to dislocations. In the following 
we will use the expression pseudomoment stress tensor for H^. 

It is known that nontrivial traction boundary problems in the variational formula- 
tion of the gauge theory of defects can be formulated by means of a so-called null La- 



grangian [12j, [l_7|, |43[ . When the null Lagrangian is added to the Lagrangian of elasticity, 



it does not change the 'classical' Euler-Lagrange equation (force equilibrium in elasticity 
theory) because the associated Euler-Lagrange equation, o - ?- ■ = 0, is identically satisfied. 
After minimal replacement, Uij 

W hg = d J ( 



Pij, the null Lagrangian 



) = a ijJ u i + °% u ^ — * a ii& 



(2.30) 



gives rise to a background stress tensor cr?- which can be considered as the nucleation field 
in the gauge theory of defects [17 . 
The Euler-Lagrange equations of W 



5W 



W Asl + W e 



dW „ dW 

- Ok- 



Wbg are given by 




6 Pij <) hj (l.Sijj,. 
which can be expressed with Eqs. ( 12. 5 p and (12.131) as 

Hijk,k + crij = cr°- (pseudomoment equilibrium). 



(2.31) 



(2.32) 



2 We want to mention that Nye [35( derived originally the Eq. ()2.25p if the elastic strain is zero (see 
also discussion in Li et al. Here we have used the differential geometrical definition (|2.24p of the 

contortion in order to derive l|2.25p . For non-zero elastic strain Eq. (|2.25|) was derived by deWit [38| . 
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Thus, the null Lagrangian gives a contribution to the pseudomoment equilibrium condi- 
tion. Because the pseudomoment stress is skewsymmetric in the indices j and k, the 
force equilibrium follows from Eq. f 1 2 . 3 2 j) as a continuity equation 

a ij,j = (force equilibrium). (2.33) 

We want to note that a similar equation like (12.321) has been earlier obtained by Trz§sowski 
and Slawianowski 44j, Kroner [45[ from a variational principle in nonlinear dislocation 
theory. We may rewrite Eq. (I2.32p as follows 

£kjnHin,k = O'ij ~ &ij (2.34) 

in order to see that the pseudomoment stress tensor Hi n is a kind of a stress function tensor 
of first order 46J. Eq. (12.321) may be decomposed into the symmetric and antisymmetric 
parts 

H(ij)k,k + = tfij) (2-35) 
H[ij]k,k + = (moment equilibrium). (2.36) 

Hence, Eq. (12.361) has an appearance close to Cosserat's moment of momentum equation 
and Huj]k has the meaning of a moment stress tensor. The divergence of the moment 
stress tensor gives rise to asymmetric force stresses. The existence of moment stresses 
H[ifl k demands asymmetric force stresses a^y However, due to Eq. (I2.35f) the equilib- 
rium equations of dislocation theory deviate distinctly from those of the Cosserat theory 
(see also Kroner 45|). Using Eqs. (12. 9p and (12.191) . we can rewrite the moment equilib- 
rium (12.321) in the following form 

Cl{@ik,jk — Pij,kk) + C2(Pji,kk ~ Pjk,ki + Pkj,ki ~ /3ki,kj) 

+ c 3 [5ij(Pik,ki - Pu,kk) + Pliji ~ Pij,ti] + °u' = a ij- (2.37) 
With the inverse constitutive relation for /3y 

7 ~\~ Li 7 — /x v 

Pij = — CTjj + — CTji — rSij (Jkk (2.38) 

4/i7 4/ry 2/41 + v) 

where the Poisson ration v is expressed in terms of the Lame coefficients 

A X = ^- (2.39) 



2(A + /i)' l-2u 

we are able to rewrite completely the pseudomoment equilibrium condition (I2.37P in terms 
of the force stress tensor according to 



x 27Z/ 

v ci - c 2 + 2c 3 ) ~ 2c 37 



(Sij (?U,kk - + 2c 3 ^8ij (Jkl,kl 

+ [01(7 + fi) - c 2 (7 - /j,)] (<Jik,jk - <?ij,kk) + [cx{l - n) - c 2 (7 + //)] (<Jki,kj - <?ji,kk) 
+ [2c 2 // - 03(7 + n)] a k j,ki - [2c 2 /i + 03(7 - //)] a jkM + 4/17 = Afij cx°-. (2.40) 
If we use the force equilibrium condition cr^-j = 0, Eq. (12.401) simplifies to 

(ci - c 2 + 2c 3 ) ~ 2c 37 ou,kk - <Tu,ij) - [ c i(7 + /^) - c 2 (7 - A*)] °ij,kk (2.41) 

+ [ci(7 - n) - c 2 (7 + /i)] {(7ki,kj ~ <?ji,kk) + [2c 2 /i - c 3 (7 + //)] (T fcjifci + 4/i7 a {j = 4/i7 <r°-. 

Equation (12.411) is the master equation, being the central equation for the force stress 
tensor in this gauge theory of dislocations. In the next two sections, solutions for screw 
and edge dislocations will be given. 
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2.2 Comparison to other gauge theories of dislocations 



At this point we want to discuss the relation to other translational gauge theories of 
dislocations which differ in the form of the constitutive relations of the force stress and 
pseudomoment stress. The differences between the gauge theories of dislocations are 
mainly the choice of the constitutive relations of the force and pseudomoment stresses. 



Edelen [111 . |12| . Il7j used a very simple constitutive relation for the pseudomoment stress 
tensor. The Edelen choice is Ci = 2s, c 2 = and c 3 = 0. Of course, this is not a 
general isotropic constitutive law. In addition he assumed symmetric force stresses. Also 
Valsakumar and Sahoo 16| used this choice for the constitutive relations. Malyshev 18 



and Lazar 2jJ, ^2fl discussed and used the so-called Einstein choice in three dimensions. 
It is given by: a 2 = —ax and a 3 = — ai/2. Thus, C\ = ai/2, c 2 = —ax/ 2 and c 3 = — ax- In 
the Einstein choice the pseudomoment stress tensor is symmetric in the first two indices 
Hijk = H(ij^k and is given in terms of 'gradients' of the elastic strain f3u,j). In this 
choice, Huj} kjk reduces to the Einstein tensor. Another choice C\ = ax/2, c 2 = — ax/2 
and c 3 = u/(l — v)a x was used by Lazar [23J to calculate symmetric force stresses of 
screw and edge dislocations. It can be seen in Eq. (12.211) that the Einstein-choice and 
also the latter choice violate the condition of nonnegative energy Wdisi > 0. The material 
coefficients ax, a 2 and a 3 have to be nonnegative. But in this paper we use the most 
general linear constitutive relations (12.61) and (12.191) of an isotropic material with six 
material parameters. In Table Q] we have listed the choice of the constitutive relations. 





Edelen [10] 


Einstein choice [18, 2L] 


Lazar [23] 


present paper 












V 


V 


p 


V 


V 


7 











7 


ax 
a 2 
a 3 


2s 
2s 
2s 


a i 
—ax 

ax 
2 


ax 
i — a-\ 

Ol 
2 


ax 
a 2 

«3 



Table 1: This table lists the material parameters o/, /i, v and 7 for the translational gauge 
theory of dislocations. 



2.3 The Green tensor 

In this subsection, we want to derive the three-dimensional Green tensor of the mas- 
ter equation (12.411) . First we set cr^- = LijS(x)S(y)5(z). The fundamental solution of 
Eq. (12.41 j) is defined by the equation 

~ 2c s7 (Sij cru,kk ~ <Tu,ij) ~ [01(7 + fi) - c 2 (7 - n)\<Tijjik ( 2 - 42 ) 

+ [ex (7 - /i) - c 2 (7 + ji)] (a k i,kj - crji,kk) + [2c 2/ u - 03(7 + //)] a kj ,ki 

+ ifijcTij = 4/ijLij 6(x)S(y)6(z). 

We use the notation for the three-dimensional Fourier transform (47I . \a 
m=^[f(r)] = r /(r)e +i ^dr, f(r) = [f(q)) = -Lj / f(q)r 



(cx - c 2 + 2c 3 



(2tt) j 



iqr dq. 

(2.43) 
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Applying the Fourier transform, it follows that 



[ci - c 2 + 2c 3 ) ~ 2c 3 7 q 2 vu ~ QiQj^ii) ~ [ci{l + fj) - c 2 (7 - AO] q 2 &ij (2.44) 

+ [ci(7 - /i) - c 2 (7 + /i)] {q k q^ki - q 2 &ji) + [2c 2 /i - c 3 (7 + At)] flfe*?*,- 

- 4/i7<x ij = -ApyLij. 

The force equilibrium (I2.33P reads in the q-space: cfijqj = 0. Symmetrization, antisym- 
metrization, taking the trace and the inner product by q from left of Eq. (12.441) we deduce 

[1 + £\q 2 ]a {ij) 1%) [SijQ 2 - qiQjH ~ (£\ - ^qk^iOj) = ( 2 - 45 ) 

[1 + ^<? 2 ]%] + (f 4 ~ ti)q k Vkm = L m (2.46) 
[1 + f 2 q 2 ]a u = L u (2.47) 
[1 + i\q 2 \q^ l3 = qi L tJ (2.48) 



where the characteristic length scales are defined by 

. 2 ci - c 2 ai 



2/i 2{jt 
[1 - v){c x - c 2 + 2c 3 ) (l-z/)a 2 



2/i(l + z/) 2^(l + i/) 

(/i + 7) (ci - c 2 + c 3 ) (/i + 7) (ai + a 2 ) 



13 4/^7 8/^7 

2 ci + c 2 01 + 2a 3 



f-4 



27 67 

The lengths £±, £2 and £ 3 fulfill the following relation 

,2 _ V + 7 ( B 2 , 1 + ^ B 2 



' 3 ~ 4 7 V^^l-z/^ 2 



(2.49) 
(2.50) 
(2.51) 
(2.52) 

(2.53) 



Now some remarks are in order. In the static gauge theory of dislocations we have found 
four characteristic length scales. These length scales may describe size effects produced 
by dislocations. They are important in the plastic zone namely the dislocation core. 
All four length scales depend on the six material parameters of an isotropic material. 
Thus, they describe the material property on which the influence of the pseudomoment 
stresses depends strongly. If the ratio of the smallest dimension of the body to the length 
scales is large, the effects of pseudomoment stresses are negligible. However, when this 
ratio is not large, pseudomoment stresses may produce effects of appreciable magnitude. 
This is the case in the dislocation core region. It is interesting to note that the internal 
length ly depends on /x and has a similar form as the internal length in the couple stress 



theory [49j, [50J and in gradient elasticity [30] - It can be seen that £2 depends on the 
Poisson number v. This length £2 is the characteristic length of dilatation. Such a 
length appears also in gradient elasticity j3o|. Because the characteristic lengths £ 3 and £4 
depend on 7 they look like the two characteristic lengths of micropolar elasticity namely 
the characteristic lengths for bending and for torsion (see e.g. 5lj]). Thus, it depends on 



the physical problem which length scale is of importance. It is clear that for 7 — >• 00 just 
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the bending length, which is the characteristic length in the theory of couple stresses, 
survives 



lim t\ 

7^00 



ci - c 2 + c 3 a\ + a 2 



4/x 



8/i 



lim . 

7—* 00 



0. 



(2.54) 



On the other hand, if 7 = 0, then £3 — > 00 and £4 — > 00. For the Einstein choice a 2 = —ax 
and a 3 = — ai/2 the lengths convert to 



2 _ £1 _ 
1 /j, 2ji' 



'1 - i/)ai 



/i(l + z/) 2fi(l + u) 



0. 



0. 



(2.55) 



Thus, only two lengths survive. These length scales £\ and £\ agree with A4~ 2 and — A/"~ 2 
used by Malyshev [18j . In addition, four length scales survive in the Edelen choice together 
with asymmetric force stresses 



[I - iAs 



'2 — , .A ' C 3 



(/^ + 7)g 
2/i 7 ' 



7 



(2.56) 



The length scales £ 2 and £\ coincide with M 2 and iV 2 introduced by Kadic and Edelen 
Tof. Edelen 52]. If 7 = 0, then ^ 3 — > 00 and £4 — > 00. 
From Eqs. (l2T45l) -( |2T48ll we obtain 



0"H = 



1 + £ 2 q 2 

^2 



1 + 



and 



L 



l + £ 2 q 2 2 



+ o(f-3)- 



1 + ^ 2 g 2 



«2 - $ 



[<% 2 ~ gjgj]^i 

a+^ 2 )(i+^ 2 : 

QkLk[iQj] 



QkLk(iqj) 



l + £ 2 q 2 )(l + £ 2 q 2 



Using the inverse Fourier transformed function [47|, 148 



^(3) 



1 



Lg 2 + 



47rr 



-r/i 



\/ x 2 + y 2 + z 2 



(2.57) 
(2.58) 

(2.59) 
(2.60) 



(2.61) 



we find 



M A7cr£ 2 



^[^A-^-][-(e-^-e-^) 
8tt lr v 7 



^— Lk(idj)dk 



1 



47rr£l 



r ^ 4 + — L k [idj]dk 



47T 



Lr 



(« 



-r/li _ -r/£3 



(2.62) 
(2.63) 



With (jjj = <T(ij) + <T[ij], the fundamental solution can be written in the following form 

GijkiLki (2.64) 



a 



1.1 
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where G^ k \ is the Green tensor of Eq. (I2.42p . So, the Green tensor reads 



G 



ijkl 



8tt 



(SikSji + Su5 jk ) 



-r/ti 



i\r 



+ (S ik 5ji - 5 u 8 jk ) 



-r/li 



(6 i:l A - didj)8 kl -(e 



ir 



Using the convolution theorem, we obtain for the solution of Eq. (12.4 ip 

&ij = Gijki * CT kl . 



(2.65) 



(2.66) 



With the help of this equation three-dimensional problems can be solved e.g. dislocation 
loops. 

3 Screw dislocation 



We consider now a straight screw dislocation. We choose the dislocation line and the 
Burgers vector in z-direction: b x = b y = 0, b z = b. First of all, a screw dislocation 
corresponds to the anti-plane strain problem. Thus, we only have the non-vanishing com- 
ponents of the elastic distortion tensor f} zx and /3 zy . From the constitutive equation (12.61) 
we obtain the relations valid in anti-plane strain: 



cr 2 



V ~7 
^ + 7 



a 



/i + 7 



a 



yz- 



(3.1) 



In the case of a screw dislocation we want to solve Eq. (12.371) directly. From the force 
equilibrium condition (I2.33p . the following condition follows 



fizx,x ~t" @zy,y 0. 



(3.2) 



It looks like a gauge condition, but it is nothing else than a consequence of the force 
equilibrium condition (I2.33p . Using Eq. (13.21) . (I2.37P simplifies to 



Cl 



/i + 7 

Cl 

/i + 7" 



-A 



A 



fizx /3 zx i 
Pzy — Pzy 



1 

1 + 



C 2 



-A 



/i-7 
/i-7 



A 



@zx 
ftzy 



0° 

r'zxi 



zy 



(3.3) 
(3.4) 



Here A denotes the two-dimensional Laplacian. Because we have four equations for two 
components (3 ZX and {3 zy , we obtain a relation between c 2 and C\ as follows^! 



(-2 



/i + 7 



Cl 



(3.5) 



3 If 7 = 0, wc obtain C2 = —c\. If C2 = like in the Edelcn choice, we obtain c\ = 0, £f = 0, f3 zx = f3® x 
and (3 zy = (3® y . Thus, only the classical solution of a screw dislocation is allowed in the Edelen choice 
with symmetric force stresses. Edelen pjj used some ad hoc assumptions like (3^ z = flzx to avoid this 
problems. But in the classical theory it must be: /3® z = 0. Also the approach of Valsakumar and Sahoo 
[lq possesses such lacks. 
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or in the 'irreducible' parameters 



«3 



37 - V 
2/i 



a 1 . 



(3.6) 



Relation (I3.5P is a consequence of (13.11) . From the relation (13. 6p together with the condition 
of non-negative dislocation core energy (I2.2ip we obtain the constraint between fi and 7: 



7> 



(3.7) 



It is important to note that in the case of 7 = 0, 03 = — cti/2 violates the condition (12.211) . 
Therefore, the condition of nonnegative dislocation core energy demands asymmetric force 
stresses. For the Edelen choice (03 = ai) we obtain from (13. 5p and (13. 6p : /i = 7. Thus, 
it does not make sense to use symmetric force stresses in the Edelen model. Due to 
these reasons, the Edelen model demands asymmetric force stresses with fi = 7 for a 
screw dislocation. Therefore, with 7 = fi the correct solution of a screw dislocation of 
the Edelen model is contained in our general solution of the screw dislocation which we 
calculate in the following. 

If we substitute ( 13.51) and ( 13. 6p in ( 12.491) and ( 12.52[) . we obtain the characteristic length 
scale of the anti-plane strain problem 



Cl 



ai 
2^ 



/i + 7 

and the equations for the elastic distortion simplify to 



2// 



0° 

Hzy 



(3.8) 



(3.9) 



They are two-dimensional inhomogeneous Helmholtz equations. The inhomogeneous parts 
are given by Eq. (1A.2|) . If we substitute (1A.2I) into Eq. (13.91) . the elastic distortion is easily 
obtained as 



Pz 



b_y_ 

2tt r 2 



1 



b x 



r 

1 - -K, 



(3.10) 



where K n is the modified Bessel function of the second kind and n — 0, 1, . . . denotes the 
order of this function. Substituting the elastic distortions ( 13.101) into ( 12.91) . the torsion 
tensor or dislocation tensor is calculated as 



zxy 



a. 



2n£i 



K n [- 



(3.11) 



The torsion ( 13. lip generates the following components of the pseudomoment stress ten- 
sor (I2.19P according 



H 



xyz 



H, 



uy 



H. 



yz.r 



H 



zxy 



C2 T-zxy 
&2 T Z xy 
Cl T ZX y 



(fi - 7)6 ^ ( r 



2tt 

2vr 
+ 7)& 



K 



2tt 



(3.12) 
(3.13) 
(3.14) 
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Thus, this localized pseudomoment stress caused by a screw dislocation is of torsion-type. 
It is given in terms of only one length scale l\. The trace of the pseudomoment stresses is 



H 



- 37)& v ( r 



(3.15) 



Now we want to show how we can solve the master equation (12.411) by means of a 
suitable stress function ansatz for the asymmetric force stresses of a screw dislocation. 
With the stress function ansatz we fulfill the force equilibrium condition (12.331) . Then we 
substitute the stress function ansatz into the pseudomoment equilibrium condition (12.411) . 
In the case of a straight screw dislocation Eq. (12.411) reduces to 



[ci(7-/i) -c 2 (7 + ^)](o r : 
[d(7-/i) - c 2 (7 + //)]( 



yz,yx 



xz,xy 



Vxz,yy) ~ [Cl(7 + //) - C 2 (j ~ //)] A(J ZX + A^ja 

&yz,xx) - [ci(7 + /i) - c 2 (7 - //)] Aa zy + Apqo 



2.'/ 



' yz,xy ) 



[c 2 (7 + /i) - 01(7 - //)] Acr za; + [c 2 (7 - /i) - 01(7 + p)]Aa xz 

+ [2c 2 /i - 03(7 + //)] (cr XZtXX + a 

[c 2 (7 + //) - 01(7 - fi)]Aa zy + [c 2 (7 - /i) - Ci(<y + fi)]Aa yz 

+ [2c 2 /i - 03(7 + //)] (a XZjXy + (y yZ) y y ) + A^a. 



(3.16) 

4/i7cr° 
(3.17) 

(3.18) 



+ 4fi~f(T xz = A^a x; 

(3.19) 



4/rycr, 



We introduce the following stress function ansatz, which is suitable for asymmetric force 
stresses in the anti-plane strain problem, 



a 



■1 j 













-d y F d x F 



-^dyF 
^d x F 




(3.20) 



where F is a modified Prandtl stress function. Using (13.51) . (12.491) and (1A.7|) . we have to 
solve the two-dimensional inhomogeneous Helmholtz equation 



Its solution reads 



[l-£lA]F = F c 
(/i + 7 )6 



2tt 



In r + Kq 



With Eqs. (I3.20p and (I3.22|) the force stress is calculated as 



0" 2 



a 



ZLJ 



(/i + 7)6 y 
2tt r 2 
(/i + 7)6 x 
2^ 7 1 



-1 



r 



1 " T« 



/ r 


)] 






r y 




TJ. 





(fi - 7)6 y 
2tt r 2 
{H~l)b x 



1 



2tt 



r 



1 - a * 



/ r 


)] 






r y 




TJ. 





(3.21) 
(3.22) 

(3.23) 
(3.24) 



In cylindrical coordinates the non-vanishing components of the force stress and elastic 
distortion tensors read 



0" 2 



Pz 



0x + 7)6 1 
2tt r 
b 1 r 
In r l 



1 " 

r ( r 



cr<t,z 



27T r 



1 " 5* 



(3.25) 
(3.26) 
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Figure 1: The component (3 ZV is given in units of 6/[27r£i] . The dashed curve represents 
the classical component. 



The force stress and elastic distortion fields, calculated in the gauge theory, do not possess 
singularities. In fact, they have at r ~ 1.114£x a maximum of: o z $ ~ 0.399(/i + 7)6/[27r£i], 
a^ z ~ 0.399(/x - 7)6/[2?r4], Pz<t> ^ 0.3996/[2vr£i] (see figured]). We want to note that if 
7 = 0, we recover the force stress and pseudomoment stress, the elastic distortion and 
the dislocation tensor earlier calculated by Malyshev [l8[ and Lazar 21, 22] in the gauge 
theory of dislocations with the so-called Einstein choice. The modified Burgers vector is 
given by 



b{r) 



2tv /•27T i*r 

(3 Z(j) rd(j)= / a zz (r')r'dr' 
o Jo Jo 



(3.27) 



This modified Burgers vector differs appreciably from the constant value b in the core 
region from r = up to r ~ 6£i. Thus, it is suggestive to take r c ~ 6£i as the disloca- 
tion core radius. Outside this core region, the Burgers vector reaches its constant value. 
Accordingly, the gauge theoretical solutions approaches the classical one outside the dislo- 
cation core region. Due to the cylindrical symmetry of a straight screw dislocation around 
the dislocation line the physical situation requires only one length scale l\. 

The force between dislocations is the so-called Peach-Koehler force 53j, |5J, |55j . In the 
gauge theory of dislocations the gauge invariant formulation of the Peach-Koehler force 
has been derived by Lazar and Anastassiadis [56( . Kadic and Edelen [ll| and Edelen and 
Lagoudas [12| just derived an expression which is not gauge invariant because they used 
the canonical Eshelby stress tensor instead of the gauge invariant one. It is given by the 
convolution of the stress and the dislocation density tensor 



— &u * Tuk = e kij a H * a>i 



(3.28) 



It is obvious that (I3.28j) is similar in the form as the classical expression of the Peach- 
Koehler force. But, in (I3.28j) the stress and the dislocation density tensor are the gauge 
theoretical ones. For a screw dislocation it reads 



F v = —a zx * a. 



(3.29) 



By the convolution theorem, it is in the Fourier space 

F x (q) = F{2) [F x ] = a zy (q)a zz (q), F y (q) = T (2 ) [F v ] = -a zx (q)a zz (q). (3.30) 
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We use the following notation for the two-dimensional Fourier transform 47], |48 

1 



The Fourier transformed functions are 

2tt 



^ (2 )[hir] = - 



Q 



(2) 



2tt 



q 2 + 



and we obtain 



a 



zy\ 



(/i + 7 )6- 



iq y 



a 2 



V- 



f(q)e-^ r dq. 
(3.31) 

(3.32) 

(3.33) 
(3.34) 



The Fourier transformed Peach-Koehler force F x is 



, 2 ' 



Using the inverse Fourier transform, the Peach-Koehler force F x is calculated as 



F x = -{fi + -f)bb'd x F- 2 ] 



q A q l 



h 



(fl + 7)66' X 



2tt 



After an analogous computation we obtain for the component F y 

_ (// + 7)66' y_ 
v 2tt r 2 



-^i(-)-^o(- 



In polar coordinates, the non-vanishing component reads: 



(/x + 7)66' 1 
27T r 



r 

1 " 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



The Peach-Koehler force between two screw dislocations is a radial force. It has a maxi- 
mum of: F r ~ 0.2488(/z + 7)6&7[ 27r ^i] at r ~ 2.324£i and is zero at r = (see figure El). 
It is interesting to note that the maximum of the Peach-Koehler force calculated in the 
framework of nonlocal elasticity is higher than our gauge theoretical result (see Eringen 



571 . |58| . Lazar [59(). If one calculates the gradient of the elastic interaction energy be- 



tween two screw dislocations, then also the result (I3.38|) follows. But, if one uses the 'total' 
interaction energy including the core contribution, one does not get the Peach-Koehler 
force. Because Valsakumar and Sahoo 6(| used the 'total' interaction energy, they have 
not calculated the Peach-Koehler force in the gauge theory. Nevertheless, their result is 
in agreement with the nonlocal expression calculated by Eringen (57| . 
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Figure 2: The component F r is given in units of fibb' /[2n£i\ (7 = 0). The small dashed 
curve represents the nonlocal result and the dashed curve represents the classical compo- 
nent. 



4 Edge dislocation 

In this section we consider a straight edge dislocation. We choose the dislocation line in 
the ^-direction and the Burgers vector reads b x = b, b y = 0, b z = 0. The problem of an 
edge dislocation corresponds to plane strain. 



4.1 The general case 

Let us turn to the edge dislocation in the framework of the gauge theory of dislocations. 
In this case the master equation (I2.37P reduces to 

[(ci - c 2 + 2c 3 ) 27Z/ - 2c 3 7(l + v)] {a xxm + o yym ) - 27(01 - c 2 )Aa xx (4.1) 
+ [(ci - c 2 )(7 -y)- c 3 (7 + (j,)] (a XX}XX + (Jy X ,yx) + 4/i7 a xx = 4/ry a° xx 

[(ci - c 2 + 2c 3 ) 27Z/ - 2c 3 7(l + v)] (a xx>xx + a yy , xx ) - 27(01 - c 2 ) Aa yy (4.2) 
+ [(ci - c 2 )(7 -/m)- c 3 (7 + y)] {(T X y jXy + v yv ,yy) + 4^7 a yy = 4/ry a yy 

[01(7 - /i) - c 2 (7 + //)] (ovrb,^ - «T WXliRl .) + [2c 2y u - c 3 (7 + 

A*)] (^"xj/j^x "I - ®yy,yx) 

- [(ci - c 2 + 2c 3 ) 2 7 z/ - 2c 37 (l + 

^)] (^11,1!/ "I - ®yy,xy) (4-3) 

- [ci(7 + /i) - c 2 (7 - /i)] Acr^ + 4//7 a xy = 4//7 a° y 
[01(7 - fi) - c 2 (7 + //)] (<T Wl j,a. - o-jp.,,^) + [2c 2 /i - C 3 (7 + A*)] {.&xx,xy @yx,yy) 

- [{c x - c 2 + 2c 3 ) 2 7 i/ - 2c 37 (l + 

^)] ^Pxx,xy &yy,xy) (4-4) 

- [01(7 + /i) - c 2 (7 - /i)] Aa p + 4/i7 cr^ = 4^7 a\ 



yx 



2/iz/ 



0",, = a 



(4.5) 



We used the plane strain condition a zz = v{a xx + a yy ) because we require plane strain for 
an edge dislocation in the gauge theory. 

For the asymmetric force stress we use the stress function ansatz of Mindlin-type with 
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the two modified stress functions / and 

( d 2 f - d 2 V -d 2 f + d 2 V \ 

yy J x v x v J xx 



V 



(4.6) 



'A/ J 



This stress function ansatz (14. 6[) has been also used by Lazar and Maugin (6l| in gradient 
micropolar elasticity. If we substitute (14. 6[) into the coupled system of equations (|4.ip - 
( 14. 5p . we obtain the decoupled equations 



dy 2 



JUL 

dx 2 



If 



(ci - c 2 + c 3 )(l - I/) 



2/i 

d 2 



1 - 



dxdy 
(ci - c 2 + c 3 )(l - 1/ 



+ 



2/i 
«9 2 



ci -c 2 + c 3 )(jU + 7) 
4/ry 

■a)/-/ " 

(ci - C 2 + C 3 )(/i + 7) 



(4.7) 



dxdy 



dxdy 
(ci - c 2 + c 3 )(l - v) 



2/2 

d 2 



4^7 
A) / - /° 



(4.t 



1 



(ci - c 2 + c 3 ) (// + 7) 



d 2 



dxdy 



dx 2 - 

(ci - c 2 + c 3 )(l - 1/) 



+ 



2/7 



4/i7 
A)./-/ 



(4.9) 



f 1 - ( C 1- C 2 + C3)(// + 7) A W ^0 

V 4/ry / 



1 - 



c 3 (l - ^ ) 
2/iz/ 



A)/-/ 



0. 



The addition of Eqs. (I4T7]) and (Oil gives 

(ci - c 2 + c 3 )(l - v) 



A 



2/^ 



A)/-/ 



0. 



For the moment equilibrium condition (12.361) . we obtain 

(ci-c 2 + c 3 )(/i + 7) A ^_^ 



1 - 



4/17 



0. 



(4.10) 
(4.11) 

(4.12) 

(4.13) 



If we compare (14. lip with (14.121) . we obtain the following constraint between ci, c 2 and 



c 3 



1 - V 



[Cl - c 2 



(4.14) 



4 This relation is not fulfilled in the Einstein choice and c\ — C2 + C3 = 0. In the Edclen choice C3 = 0, 
c 2 = 0, we obtain from lj4.14jl ci = 0, li = and ^3 = 0. Thus, only / = f° and * = *° are allowed with 
the Edelen choice if the plane strain condition is fulfilled. Due to these reasons, Edelen [EH, Kadic and 
Edelen 11 1, Malyshev [18[ used a stress function ansatz not fulfilling the plane strain condition. 
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or in 'irreducible' parameters 



«2 



1 + 1/ 
1 - V 



(4.15) 



This constraint is nothing but a consequence of the relation a zz = v{a xx + a yy ) valid in 
the plane strain problem. We want to note that this constraint (14.151) has been earlier 
derived by Lazar 23J. With — 1 < ^ < 1/2, we get from (14.151) the relation 



< d2 < 3ai. 



(4.16) 



If we substitute (I4.14p and (I4.15P into (12.501) and (12. 51 j) . we obtain the two characteristic 
length scales of the plane strain problem 



2 Ci — c 2 i 
2 2/2 2/i 
, 2 (ci-c 2 )(m + 7) 



/i + 7 



d 4/x(l - j/)7 4^(1 - z/)7 
The relation between these internal lengths is 



CL\. 



V + l 
2(1-% 



(4.17) 
(4.18) 

(4.19) 



Therefore, we obtain from Eqs. (14. 71) - (14. lip two inhomogeneous Helmholtz equations 
which we have to solve 



[l-£ 2 2 A]f = f° 
[1 - t\ A]^ = ^° 



(4.20) 
(4.21) 



where the inhomogeneous parts are given by ( 1B.14I) and ( IB. 151) . The solutions of Eqs. (14.201) 
and (I4.2ip are the following modified 'Airy' stress functions 



./ 



fib 



4tt(1-i/) y 



d v {r 2 \nr + 4£ 2 2 



d x {r 2 hir + 4£ 2 3 



lnr + iToff 

V *-2 

lnr + K (^ 

\l 3 / j 



(4.22) 
(4.23) 



If we substitute the stress functions (14.221) and (14.231) into the stress function ansatz (14. 6p . 
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the following components of the force stress tensor follow 



(To 
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[X 



(y 2 + 3s") + ^% 2 - 3x 2 ) - 2y*-K 1 - - 2(y a - 3z 2 )K 2 
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_ ^) _ IIS^ _ 3x 2 ) - 2* 2 -^f 1 ( -) + 2(y> - 3x')K 2 (- 



■1' J 
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4<! A 



2^ 4£ 2 

4. 



(3x 2 -y 2 ) + 2x 2 -K 1 

£3 



+ 2(3x 2 - y 2 )K 2 



x 2 - y 2 ) --fix 2 - 3y 2 ) - 2y 2 -K 1 - + 2(x 2 - 3y 2 )K 2 - 



B 



[x 2 + 3y 2 ) + ^(x 2 -3y 2 ) 



2x 2 ^-K l 



2{x 2 -3y 2 )K 2 



a 



4<M 



r 

- B 



4> 



, x- 



y 



^(x 2 - 3y 2 ) - 2y 2 -K 1 [-)+ 2{x 2 - 3y 2 )K 2 ( - 



(x 2 -y 2 )-%x 2 -3y 2 )-2y 2 ^K/' 
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■3' J 



with 
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1 ~tMt 2 



A :-- 



fib 



2tt(1 - i/l 
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^76 



The trace of the stress tensor is 



O'kk 



-2(l + u)A^- 



vr(/i + 7)' 



The skew-symmetric part of the force stress tensor reads 



a 



/ryfe a; 



7r(/i + 7) r 2 



5*G 



(4.24) 



(4.25) 



(4.26) 



(4.27) 
(4.28) 

(4.29) 
(4.30) 

(4.31) 



Now we want to discuss some details of the core modification of the force stress fields (14.241) - 
(14.28p . The spatial distributions of the force stresses near the dislocation line are presented 
in figure [31 The force stress fields have no artificial singularities at the core and the max- 
imum stress occurs at a short distance away from the dislocation line (see figures H] and 
[5]). They are zero at r = 0. It can be seen that the stresses have the following extreme 
values: \a xx (0,y)\ ~ 0.546A/4 + 0.2605/4 at \y\ ~ (0.996£ 2 + 1.494£ 3 )/2, |<r w (0,y)| ~ 
0.260^4/4-0.2605/4 at \y\ ~ (1.4944 + 1.494£ 3 )/2, \cr xy {x, 0)| ~ 0.260A/4 + 0.5465/4 
at \x\ ~ (1.4944 + 0.9964)/2, \a yx (x,0)\ ~ 0.260A/4 - 0.2605/4 at \x\ ~ (1.4944 + 
1.4944) A and \a zz (0,y)\ ~ 0. 399,4 at \y\ ~ 1.1144- Thus, the characteristic internal 
lengths £ 2 and £3 determine the position and the magnitude of the stress maxima. For 
7 > the stresses a xx and a xy are bigger than in the case 7 = and a yy and a yx are 
smaller than in the case 7 = (see figure [5]). 
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Using Eq. (j2.38p . we find for the elastic distortion 
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2{y 2 - 3x 2 )ir 2 
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i|(3x 2 -y 2 ) + 2x 2 fir/ r 



2(y 2 - 3x 2 )K 2 
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The dilatation reads 
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In addition the elastic rotation is 
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In Eqs. (14.361) and (14.371) it can be seen that l 2 and £3 are the characteristic lengths for 
the elastic dilatation and elastic rotation, respectively. 

Using the elastic distortion (I4.32l) -( 14.35l) in terms of the stress functions / and 
we obtain for the torsion (12. 9ft or the dislocation density of a straight edge dislocation 
according 
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Simple differentiation gives the non-vanishing expressions 
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a 



b xy 
27r r 2 



l 3 



-3' J 



(4.40) 
(4.41) 



So far, it is surprising that the component (14.41 j) . which is usually the dislocation density 
of an edge dislocation with Burgers vector b y , is non-zero. The components ( 14.401) and 
( I4.4ip are necessary to fulfill the pseudomomentum equilibrium condition f!2.32[) . Also 
we would like to note that these non-vanishing components of the torsion tensor do not 
have cylindrical symmetry due to the .KVterms (see figure [6]). Since an edge dislocation is 
lacking cylindrical symmetry around the dislocation line two length scales, £2 and £3, are 
needed for a proper model. Fig. Eh shows this asymmetry in the dislocation core region 
of an edge dislocation. 

With Eqs. (I2.19P and (12.22p the localized pseudomoment stresses of bending type are 
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Figure 4: Force stress of an edge dislocation: (a) a xx , (b) a xy , (c) <j yy , (d) a yx are given 
in units of A with v = 0.3 and 7 = /i/2. 
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(4.43) 

\ 

(4.44) 
(4.45) 



Therefore, plane strain bending is given in terms of the two length scales t 2 and £3. If 
we use the components (I4.40p and (14.411) of the dislocation density, the Burgers vector is 
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Figure 5: The force stresses components near the dislocation line: (a) a xx (0,y), (b) 
&xy(x, 0), (c) a yy (0,y), (d) <Jy X (x, 0) are given in units of A with z/ = 0.3 and 7 = yu/2. 
The dashed curves represent the stresses in asymmetric elasticity and the small dashed 
curves the symmetric force stresses (7 = 0). 
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Figure 6: Torsion contours of an edge dislocation near the dislocation line: (a) a xz , (b) 
a yz with v = 0.3 and 7 = ji/2 (in units of 6/[47r]). 
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Figure 7: The modified Burgers vector of an edge dislocation b(r)/b for £2 = £3 (small 
dashed) and for v = 0.3 and 7 = /i/2 (solid). 



calculated as 
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(4.46) 
(4.47) 



Thus, it can be seen that the torsion (I4.4ip does not contribute to the Burgers vector. 
Only the i^o-terms in (14.401) give a contribution to the Burgers vector (14 .46 \\ . The plot is 
given in figure H 



4.2 The case £ 2 = £3 

Now we turn to a special case of the gauge theoretical edge dislocation we have found. We 
obtained two internal lengths £2 and £3. If we require that the dislocation density ( 14.411) 
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has to be zero, the internal lengths must fulfill £2 = £3- Therefore, we may reduce them 
to only one internal length. If the internal lengths £2 and £3 are equal, we find 
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(4.48) 



Thus, in this case the parameter 7 is expressed in terms of /1 and v. Eventually, the 
non- vanishing components of the stress tensor are 
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Thus, a yx = and o yy = 0. The non-vanishing components of the elastic distortion are 
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Because the components (3 yx and (3 yy are compatible, they can be expressed in terms of a 
displacement field u y according to 
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Let us compare this type of edge dislocation with the edge dislocation originally introduced 
by Volterra (see also 66]). The displacement field of a Volterra edge dislocation reads 
u x = 6/ (27r) arctan|//x and u y = 6/(27r)lnr. But the corresponding force stress of this 
displacement field calculated in the theory of (symmetric) elasticity does not fulfill the 
force equilibrium condition because it produces line forces at the dislocation line and the 
dilatation and the hydrostatic pressure are zero. These drawbacks are not acceptable 
from the physical point of view. In the limit £2^0, our displacement field u y has a pre- 
factor v/(l — v) instead of 1. But the modified pre-factor is necessary to satisfy the force 
equilibrium condition with asymmetric force stresses and to give the correct dilatation. 
Thus, we have found the gauge theoretical version of an edge dislocation of Volterra type 
which is able to remove the original drawbacks of the edge dislocation of Volterra type. 
The limiting value v — 1/2 defines incompressibility and we obtain the pre-factor 1. In 
this limit A and 7 tend to infinity. Thus, in the limit of incompressibility the Volterra 
dislocation is also valid. 

The non-vanishing component of the dislocation density tensor is given by 



a., 



(4.57) 
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Thus, this dislocation density has cylindrical symmetry and has the same form as a screw 
dislocation (13.111) . The localized pseudomoment stress reads 



fib 
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vfib (r 

rizx 7z r iVo I y 
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(4.58) 



It is the pseudomoment stress of an edge dislocation with cylindrical symmetry (see 
also [23]). The corresponding Burgers vector is 



b(r) 



r ( r 
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4.3 The case 7 = 0: symmetric force stresses 

Now we want to revert the case of symmetric force stresses. If we set 7 = (symmetric 
force stress), we obtain from Eqs. (1331) . (l2T4^j) . (Oil) . (T23Uj) and (12311) : 
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(4.60) 



It is obvious that in (14.601) violates the condition of positive dislocation core energy 
( I2.2ip . Thus, the condition of nonnegative dislocation core energy requires asymmetric 
force stresses. The condition (I4.60p has only one independent internal length and it is 
the choice introduced by Lazar |23| . With B = we recover from Eqs. (14.241) — (14.28 1) the 



force stress earlier calculated by |23J and [63J, [64J in the framework of a gauge theory 



with symmetric force stresses and strain gradient elasticity, respectively. Also the distor- 
tions (14.321) — (14.35 1) become symmetric with B = and the elastic rotation (14.371) becomes 
zero. 



4.4 The case 7 — > 00 

If we set 7 -> 00 , then we obtain from Eqs. (I3T5D . (l2T49|) . (14141) . (l2T50|) and (123T|) : 
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This case does not violate the positive defmiteness of the dislocation core energy (12.211) . 
So we obtain the relation between l 2 and £3 
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The skewsymmetric stress ( 14 . 3 1|) converts to 
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(4.62) 



(4.63) 



and the skewsymmetric distortion (I4.37I) vanishes (3\ xy \ = 0. All other quantities can be 
easily calculated with 04. 6 II) . Because the elastic rotation is zero in the present case the 
corresponding energy converts to 



lim W(f3, curl/3) = W(e, curie) = W(e, curie 

7— >oo 
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(4.64) 
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where e and e p denote the elastic and plastic strains. Such a theory looks like a theory 
of strain gradient plasticity (see, e.g., Gurtin and Anand |65|.) In such a theory the trace 
of the dislocation density tensor is zero 

a ij = ~ e m e m = £jkie%j = 0. (4.65) 
Thus, it is a theory with two constraints 

7 -> oo : %] = 

a 3 — > oo : ctjj = 0. (4.66) 

5 Conclusion 

In this paper, we have investigated the (static) translational gauge theory of disloca- 
tions. We have used the most general linear and isotropic constitutive relations between 
the force stress and elastic distortion tensors and the pseudomoment stress and disloca- 
tion density (torsion) tensors. Thus, the linear theory possesses six material coefficients 
and the force and pseudomoment stress tensors are asymmetric. We have found four 
characteristic length scales for the gauge theory of dislocation given in terms of the six 
material coefficients to account for size effects of small-scale problems. We have derived 
the conditions of positive energy for the six material coefficients. We observed that the 
so-called Einstein choice of the three coefficients of the pseudomoment stress tensor vio- 
lates the positive definiteness of the dislocation core energy. This fact demonstrates that 
the three-dimensional Einsteinian gauge approach is not a realistic model for dislocations. 
We have observed that in the gauge theory of dislocations moment stress is the response 
to contortion and pseudomoment stress is the response to pseudomoment stress. Using 
Fourier transform, we have calculated the three-dimensional Green tensor of our so-called 
master equation of the force stress tensor. In the case of generalized plane strain the four 
characteristic lengths reduce to only two lengths and the number of independent material 
coefficients simplifies to four. Later, we have solved the anti-plane strain problem of a 
screw dislocation and the plane strain problem of an edge dislocation in the framework of 
gauge theory. In this turn, we have found new solutions for all physical state quantities 
of a screw dislocation and an edge dislocation. 
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A Screw dislocation in asymmetric elasticity 

In asymmetric elasticity the force stress tensor is asymmetric and fulfills the force equi- 
librium condition 

<, = 0, 4^4. (A.l) 
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The elastic distortion of a Volterra screw dislocation is given by 66 
0° -- 
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27T r 2 ' 



3° — — 



2 2,2 

r — x + y . 



(A.2) 



Thus, the present problem is of anti-plane strain type. The non-vanishing component of 
the dislocation density of a screw dislocation reads 



_ aO _ nO 

^zz r^zy,x Pzx,y 

It is the incompatibility condition. In terms of force stresses it reads 
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For a straight screw dislocation it has the form 

a° zz = b5(x)5(y). 
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In terms of the Prandtl stress function F°, the incompatibility condition is given by 



AF° = (n + j)b5(x)6{y). 



(A.6) 



The Prandtl stress function is nothing but the Green function of the two-dimensional 
Laplace equation. It reads [39 
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B Edge dislocation in asymmetric elasticity 

In the case of straight edge dislocations the equations of incompatibility take the form 
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Expressing the elastic distortions in terms of force stresses and using a^j = and a 

v WL + cx ™) ; we obtain 
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Because we deal with asymmetric force stresses we use the stress function ansatz given 



by Mindlin for couple-stress theory [50 
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where f° and \I/ are stress functions of second order. The stress function ansatz (IB. 71) 
is the generalization of the stress function ansatz with the Airy stress function f° for 
symmetric stresses. If \l/ is zero, (1B.7|) reduces to the usual expression for the stresses 
in terms of the Airy stress function f°. Equations ( IB. 51) and ( IB. 61) are reduced to the 
following inhomogeneous bi-harmonic equations 
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Because we want to consider a straight edge dislocation with the Burgers vector b = b x , 
the dislocation density tensor has the form 



So we obtain 
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Since the two-dimensional Green function of the bi-harmonic equation is 
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the solutions of (IB. Ill) and (IB. 121) are the following Airy stress functions 39 
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(IB. 14[) is the well-known Airy stress function for an edge dislocation with Burgers vector 
b x and (1B.15j) looks like an Airy stress function for an edge dislocation with Burgers vector 
by with a different pre-factor. 
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